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This paper considers the propagation of scalar (acoustic) waves from 
a single-frequency point source imbedded in a medium with random re- 
fractive index, in contrast with the usual plane-wave case in which the 
source is far removed from the medium. With the index being a statistically 
homogeneous and isotropic function of position, but not a function of time, 
the average complex field u (r) = (u(r)> and the spatial covariance 
<u,(r)u*(p)> of the fluctuation field Ui(r) = u(r) — u (r) are calculated. 
Beyond a few correlation lengths from the source, the average field can be 
approximated by a spherical wave with the same complex wavenumber 
found in the plane-wave case. A near-source wave number is also obtained. 
Under an improved far-field condition, the spatial covariance is reduced 
to spectral integration formulas for both transverse and longitudinal separa- 
tion of the receiving points. These formulas reveal that correlation lengths 
are much longer in the point-source case than in the plane-wave case, even 
though the relative variances are the same. We illustrate this result with 
plots for an exponential index spectrum and for a constant spectrum. 

I. INTRODUCTION 

For analysis of a detection or communication system which proc- 
esses signals from an array of sensors, a convenient postulate is that 
the signal field in the vicinity of the array is a plane wave (or perhaps 
a finite collection of plane waves in the multipath case) . Under such a 
postulate, coherent addition of the sensor outputs can yield array 
gain and directivity in the presence of ambient noise. However, there 
is always some disparity between the predicted performance and the 
performance realized in practice. In part, the disparity can be at- 
tributed to shortcomings in the signal model, the field not being a 
time-invariant plane wave in the vicinity of the array. The output of a 
single sensor may not be constant in time but instead is apt to 
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fade. Moreover, the outputs of different sensors do not fade "in step"; 
that is, after the array is steered, the signals do not fade with the 
unity correlation predicted by a fading-plane-wave model. Instead, 
the signals fade with correlation less than unity. The origin of these 
fading phenomena is the subject of this paper. 

A simplified model of fading is considered within the framework 
of the following assumptions. For a short period of time, the trans- 
mission properties of a propagation path are constant; then they 
undergo small deviations to attain another constant configuration for 
the next short period of time. These short-term deviations are relative 
to some nominal or average configuration, as opposed to representing 
a slow gross trend of the overall path properties. Such short-term 
deviations are modeled here by the effects of random fluctuations of 
the index of refraction, which could be associated in the under- 
water acoustic case, for example, with the temperature microstructure, 
turbulence, and circulatory motion of water masses. Deviations of path 
properties associated with fluctuations of a surface of reflection are 
not incorporated into the model. Thus, the model is most appropriate 
for short-term deviations of the properties of a pure-refracted path. 

In the specific situation analyzed below, the acoustic source is a 
single-frequency point source suspended far from any boundaries. If 
the refractive index were nonrandom and not position dependent, 
the acoustic field would be the usual spherical wave. Instead, the re- 
fractive index is a random function of position, but not of time. The 
average value of the index is not position dependent, so that the aver- 
age line-of-sight ray path is straight rather than bent. The spatial 
covariance of the index is a function of the magnitude of the position- 
difference vector (the index is second-order homogeneous and iso- 
tropic) . The problem is to find the average and spatial covariance of 
the acoustic field. 

Much of the literature (for example, Refs. 1-3 and most of Ref. 4) 
is concerned not with the above spherical-wave problem but with a 
situation in which plane waves impinge upon a half-space with 
random refractive index. One essential difference is that the spherical- 
wave source is imbedded in the random medium whereas the plane- 
wave source is far removed from the random medium. Regardless of 
how large the distance from the spherical-wave source to an observa- 
tion point becomes, this difference of configuration is preserved. f 



t The configurations are called the "radio link problem" (spherical) and the 
"radio star problem" (planar) in Ref. 9. 
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Some aspects of the spherical-wave case have been treated with the 
Rytov method (Refs. 4^5) and other techniques (Refs. 6-8) . 

This analysis treats the spherical-wave problem with a version of 
perturbation theory previously applied to the plane-wave problem. 2 ' 3 
For distances greater than a few correlation lengths, the average 
field can be approximated by a spherical wave with the same complex- 
valued wave number previously derived in the plane-wave case. 2 A 
near-source wave number is also obtained. On the other hand, it is 
found that the covariance function of the fluctuation field exhibits 
much larger correlation lengths in the spherical-wave case than in 
the plane-wave case. 3,4 This conclusion follows from simple integra- 
tion formulas for the covariances and is illustrated by plots of the 
covariance for special cases. 

II. PERTURBATION THEORY 

We consider the propagation of acoustic waves in a random time- 
invariant medium for the case of a monochromatic omnidirectional 
source. Our model is the Helmholtz equation: 

[V a + (1 + n(r))*K]u(r) = -8(r) (1) 

where V 2 is the Laplacian V- V, n(r) is the random deviation of the 
index of refraction which is a function of position r, k 2 = w 2 /c 2 , c is the 
sound velocity for a homogeneous medium if n were everywhere zero, 
co is the angular frequency of the source, u(r) is the complex amplitude 
(for example, the displacement potential) , and 5 is the Dirac delta func- 
tion. The time dependence exp (—iut) has been suppressed. We assume 
the source is suspended far from any boundaries; that is, we consider the 
medium to be unbounded. 

Our interest is in both the mean field (u) = u c (coherent field) and 
the fluctuation field u — u c = tm (incoherent field), where ( ) de- 
notes expectation. 

We develop a pair of equations for u c and Ui as follows. Consider 

[L + eL, + e 2 L 2 ]u = / (2) 

where L is a linear deterministic operator, L x and L 2 are linear 
stochastic operators, e is a size parameter, and / is a deterministic 
forcing function. With /* in (1) replaced by £ /x, the correspondence of 
(1) and (2) is evident. We put u = u c + v* into (2) and operate with 
( ) to obtain 

[L + €<L,> + e 2 (L 2 )K = / - e(L lWl > - e 2 <L 2 M,). (3) 
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We then subtract (3) from (2) in which u = u c + u t to obtain 

[Lui + € (L,m, - {LyUi)) + e 2 (L 2 w, - <L 2 m,»] 

= -e(L, - (£,))«. - e 2 (L 2 - (L 2 ))u e . (4) 

Equation (3) shows the source / of the mean field is countered by the 
sink e(L!Ui) + e 2 (L 2 Ui) describing the effects of scattering into the 
fluctuation field. Equation (4) is not written to exhibit true sources of 
Ui as much as to exhibit a zero-mean forcing function and zero-mean 
terms on the left side. These equations are generalizations of those 
derived by Keller [Ref. 2, p. 166, equations (12) through (13)] for 
other purposes. 

Solution of (3) and (4) can proceed with perturbation theory for 
the case of small e. Relative to e -» 0, equation (3) exhibits u c = 0(1) 
and (4) exhibits u t = 0(e). Accordingly, (3)-(4) can be rewritten 

[L + e(Ly) + e 2 (L 2 )K = / - e(L lUi ) + 0(e 3 ) (5) 

Lu< = -e(L, - (L^K + 0(e 2 ). (6) 

These equations can be partially uncoupled by operating with L _1 on 
(6) and substituting into (5) to obtain 

[L + e(L.) + e 2 (L 2 ))u e = / + e^LiL^LJu. 

-e^L-^L^u. + 0(e 3 ) (7) 

Ui = -eL" 1 ^ - (LOK + 0(e 2 ). (8) 

Equation (7) for the mean field u c is the result obtained by Keller 
(Ref. 2, p. 148. equation (10) ), who used a successive-substitution 
solution of (2) in conjunction with a crucial, and at-first-glance 
mysterious, replacement of L" 1 / by (u). Equation (1.8) is a version of 
Keller's equation (31) on p. 169 of Ref. 2. Thus, we have shown that 
these equations arise quite naturally from the pair (3) and (4). 

We now specialize (7) and (8) to the case of the Helmholtz 
equation (1). Here, 

L = V 2 + kl L x = 2 M (r)/c 2 L 2 = n\r)k 2 . (9) 

We assume (/*(r)) = 0; that is to say, we neglect any systematic 
dependence of refractive index upon position (the average profile). 
We have 

L-'s=-I^ r h ^T l) ° (r ' )dr '' (10) 
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where the integral is over all space. The inverse L _1 is an integral 
operator with kernel corresponding to the Green's function 

Ofr,rO-- M P [ ^ |r ~/ |] - (ID 

4x | r — r I 

Thus, the pair (7) and (8) is specialized to 
[V 2 + k» + eV(r)»K(r) 

- -8(r) - 4e 2 kl f ^P^lr-r'l] < M ( r)M ( r ')X(r') dr' + 0(e 3 ) 
J 47r | r — r | 

(12) 

u,(r) = 2elc] [ g ' [ *' ' r ~ ? f l] p(r'K(r') ^ + 0(e 2 ). (13) 
w J 47r | r — r | 

III. THE AVERAGE FIELD 

We now develop an approximation of the solution of (12) for the 
average field u c . It is assumed that the refractive index is statistically 
homogeneous and isotropic. The index covariance function is 

r(|r-r'|) - (mWm(O). (14) 

Equation (12) becomes 

{V 2 + /c 2 [l + 6 2 r(0)]}w e (r) 

. _ B(r) - 4e 2 /c 4 [ g^feJ-gJ] r( | p | )tte(r + P ) rfp + 0(e 3 ). 

(15) 
We assume an approximation of u c (r) of the form 

exp [ik 1 r |] ^ 

4t | r | 

where fc is a constant wave number to be determined (k j£ k ). It will be 
found that (16) is not a global solution, because a constant k cannot 
exist. Nevertheless, (16) can serve as a useful local approximation of 
the solution, with k interpreted as a weak and slowly varying function 
of | r |. 

If (16) were the solution, then u c would satisfy 

[V 2 + k 2 ]u c (r) = -5(r). (17) 

Then (15) and (17) yield 
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\kl[l + e 2 r(0)] - k 2 }u c (r) 

= -4eX f QXp[ik ° ' P ' ] r(| p \)u c (r + P )d P + 0(e 3 ). (18) 

The volume integral in (18) can be evaluated by an integration over 
the surface of a sphere with radius R followed by a radial integration 
from R = to R = oo. For the surface integration, we need only 
observe 

( M(r+ ) dS 2 

J S:\p\-R ' 4x2?" 

« fl (r) 8 -^, 0<R<\r\ 

kR (19) 

. . . . | | exp [ik(R — | r |)] D , , 

u e (r) sin fc J r I — , „ — ' — — , R > \ r \ 

where dS is a differential of area on the sphere S = {p: \p\ = R}. 
This mean value theorem follows from (16) and (17) ; see Appendix 
A. Then (19) inserted into (18) yields 

(fc a -tf[l + e a r(0)]}M c (r) 

aJjc 4 T r lrl 

= ^7-^ u c (r) / exp (ik R)T(R) sin kR dR 
Ic _J 

+ J exp (ik R)T(R) sin k \ r \ exp [ik(R - \ r |)] dR + 0(e 8 ). 

(20) 

If (16) were an exact global solution, then u c (r) could be cancelled 
in (20) ; the result would be a relation for the supposedly constant 
wave number k. But the integrals in (20) suggest that the relation is 
| r [-dependent, which is a contradiction. Nevertheless, (16) will 
serve as a local approximation of u c {r) in regions in which k is vir- 
tually constant. 

The following manipulations are made upon the integrals in (20). 
We run the first integral from to °o and correct for its contribution 
from | r | to oo by another term in the second integral. We then change 
the variable of integration of the resultant second integral. Then (20) 
becomes 



f <[f. 



k 2 = k 2 (l + e 2 r(0)) + ^p | / exp (ik K)V(R) sin I.-R dh' 
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- exp [i(k -k) \r\] f exp (tfe„K)r(| r\+R) sin kR dR \ + 0(e 8 ) 

(21) 

The | r | -dependence is now confined to the second integral. The 
large- 1 r | case occurs when we can assume this integral to be negli- 
gible, namely 

|exp [i(k B - k) | r |]| r(| r | + R) « T(R), R t [0, R ], (22) 

where we assume the first integration can run from to R with 
little error. The condition (22) shows that | r | must be much larger 
than a correlation distance; moreover, (22) shows that the increasing 
function exp[(Im k) \r\] must be taken into account. 
Thus, for large I r I, the wave number k satisfies 



k 2 « k 2 [l + e 2 r(0)] 

4e 2 k* 



+ 



k 



r exp (jhoRWR) sin kR dR + 0( € 3 ) . (23) 

•'n 



This is the relation found by Keller [Ref. 2, p. 151, equation (14) ] for 
the plane-wave problem. As expected, the spherical wave solution far 
from the source has the same wave number as the plane-wave solu- 
tion. 

The small- 1 r \ case occurs when the integrals in (21) nearly cancel 
one another ; the wave number k is given by 

k 2 f* fc»[l + e 2 r(0)] + 0(e 3 ) (24a) 

or 

k « fc [l + i e 2 r(0)] + 0(e 3 ). (24b) 

Whereas (24) yields the small- 1 r \ values of k directly, notice that 
(23) determines the large- 1 r \ values of k in an implicit fashion. 
However, an explicit approximation of the large- 1 r | value of k can 
be obtained. Notice that (23) could be solved by successive sub- 
stitutions, the first step employing either k in the integral below or 
employing (24) as follows 



4e 2 fc 

IT Jo exp ^ 



f exp {ik R)T(R) sin kR dR 

Jo 

fe„[i + ?i 2 r(0)] Jf exp (t ' fc ^ )r(i2) s ' m {k ° [1 + ^ 2 no)]R} dR 
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tt 4eX f exp (ik R)T(R) {sin (Jk.fi) cos fc o [^ 2 r(0)fl] 

+ cos (kji) sin k o [%e 2 T(0)R]} dR 
« 4e 2 & 3 / exp (ikJtyriR) sin k R dR. (25) 

Jo 

Since terms have been discarded consistently insofar as powers of c 
are concerned, approximations (23) and (25) yield 

k 2 fiH k 2 [l + e 2 r(0)] 

+ 4e 2 & 3 [ exp (ikJZWR) sin k„R dR + 0(e 3 ). (26) 

Jo 

From (26) , it follows that 
k « fc„[l + i e 2 r(0)] 

+ 2e 2 fc 2 f exp (ikJIWR) sin kji dR + 0(e 3 ), (27) 

Jo 

or equivalently 

Re k « /c„[l + £e 2 r(0)] + e 2 k 2 „ f° T(R) sin 2/c tf dR + 0(e 3 ), (28) 
and 

it 00 

Im k « e'/c 2 / (1 - cos 2k R)T(R) dR -f 0(e 3 ). (29) 

If r has a correlation length L and if k L X> 1 (a large-scale condition 
not yet imposed) , then 



Im k ^ e 2 # 



/" T(R)dR. (30) 

«/o 



Also, accuracy of the approximation (25) requires the bracketed fac- 
tor in the integrand to be equivalent to sin k R ; this holds when 

k o e 2 T(0)L o « 1. (31) 

But 

f T(R) dR ~ r(0)L o , 

Jo 

and (31) is equivalent to 
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Im k « k . (32) 

When (32) is not met, neither (28) nor (30) can be expected to be 
a good approximation. Also, for the successive-substitution procedure 
to yield a good approximation at this first step, it appears sufficient 
that the first step value (28) be well approximated by the initial 
value (24) ; equivalently, 

r(0) » k„ f T(R) sin 2k R dR, (33) 

which is a restriction on the large-wave number value of an integral 
which resembles the spectrum of T. 

The approximation (16) for the average field u c {r) together with 
(23) and (24) for the large- 1 r | and small- 1 r | values of k comprise 
the principal results of this section. The further approximations (28) 
through (30) for the large- 1 r | case are more useful than (23), but 
conditions (31) through (33) must be met. When (28) through (30) 
are compared with the small- 1 r | approximation (24), it can be seen 
that the spherical wave (16) develops attenuation and a change in 
phase velocity as \r\ increases. The transition from small-|r| to 
large- 1 r | behavior occurs when (22) begins to hold, namely, when 
the second integral in (21) begins to become negligible. The order of 
magnitude of this transitional value of | r | is a few correlation lengths. 

IV. COVARIANCE OF THE FLUCTUATION FIELD 

The previous section provides a solution of (7) or (25) for the av- 
erage field u c (r) which now can be used in (8) or (13) to yield the 
fluctuation field w f (r). Thus, (11) and (13) yield 

U<(r) = -2efc 2 a f G(r, r'Mr')u c (r') dr' + 0(e 2 ), (34) 

where u e (r') is given by (16) in which k is a weak function of \r' \. 
The spatial covariance function (w,(r)w*(p)) is now computed for the 
case in which the medium is statistically homogeneous and isotropic. 
Equations (34) and (14) yield 



<t*.(r)tt*,(p)> = 4e 2 k 2 jj G(r, r')G*( P> p') 



•r(| r' - p' \)u c {r')u*(p') dr' dp' + 0(e 3 ). (35) 
It is convenient to change to the following variables of integration 
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(with unity Jacobian) : 



where 



r' + p' 



r = y + o 



x = r' - p' , 



p' = y- 



(36) 



(37) 



Moreover, it is convenient to evaluate the fields at the following 
points: 



r-, + | 



P = v 



t 

2 ' 



where, by definition, 



V = 



r + p 



£ - r - p. 



(38) 



(39) 



The relation of the positions (36) -(39) is shown in Fig. 1. The 
covariance of the fluctuation field w< is thus 

U<(v + |)w*(t7 - ^ 

= tfk\ JJ g( v + 1 , y + §)<?*(u - I , y - |)r(| * |) 

•«.(» + f)u*(fl - l) dx dy + 0(e 3 ). (40) 

In words, the second-moment of the fields at observation center -q 
with observation position-difference vector £ comprises the integrated 
effect of scattering of the average field by the refractive index at scat- 
tering center y with scattering position-difference vector x. 

We now approximate the integrand of (40). Although the approxi- 




Fig. 1 — Scattering points f, p' and receiver points r, p. 
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raations are not valid over all space, they are valid for a region which 
can account for the major contribution to (40) in the case to be de- 
scribed later. The first approximation involves 



G(r, r')G*(p, p') = 



exp [ik (\ r — r' \ — \ p — p' |)] 

It \2 I 'II 'I 

(4tt) \r —r \-\ p - p | 



(41) 



But 



and 



r-r'| = h- 2/ + K*-*)| 

p — p' I = \v — y — !(£ — x)\ 

i „. i K^ - g)-(^ - y) 

= \ v — v \ i —\ 



+ 



(42a) 



(42b) 



The above expansions in powers of (£ — x) are appropriate for a large 
vector 77 — y as perturbed by the small vectors ±£(£ — a;). Our approxi- 
mation of (41) is 



' ;[i > + 2' y + lr v~2' y 



°^fe^ 



V I 



(4 7 



v - y 



exp [ik,(y)-(Z - x)] 
(4tt) 2 I v - y I 2 



(43) 
(44) 



where the relation 



k.(y) = K 



(45) 



defines a scattering wavevector . 

The second approximation involves replacing the coherent-field factor 
u e (y + x/2)u*(y — x/2) by a function that locally represents the 
fields as plane waves. Thus, it follows from (16) that 



« \y + %y*c[y -§) = 



exp 



ill 



2 



- k* 



(4*) ! 



i x 

y + 2 



y: 



(46) 
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where the wavenumber k is a weak function of position. But 







X 


y 


+ 


2 

X 


V 




2 



i , i y 



= I y 



■x + 



2\y 



leads to the approximation 



where 



My + 'o) u *c\y ~ g 



u c (y) 



k(y) = (Re k) 



defines an incident wavevector, and 

I ,. /.a ia _ exp(-2Imfc \y |) 

Collecting these approximations into (40) yields 



(47a) 
(47b) 

(48) 

(49) 

(50) 



uAv + Mv-l 



4e 2 /c 



o J dy 



exp [ik„(y) •£] exp [— 2 Im k | y |] 



(4t) I v - y I I y I 
•J daj T(| jc |; exp {i[k(y) - k.(y)]-x}. 



(51) 



This is the central result of this section. Equation (51) has the 
physical interpretation of a volume distribution of sources. The 
source at y generates a plane wave at the receiver with correlation 
exp[ik s (y) •£]. The strength of this wave is proportional to | -q — y \~ 2 
| y I -2 and to the value of the spectrum 



S(\ k |) = J dx r(| x |) exp {ik-x) 



(52) 



as evaluated at the local wavevector k(y) — k s {y). At this wave 
vector, the spectrum is a measure of the amplitude of those com- 
ponents of refractive index with the orientation and the periodicity 
required for constructive interference (Bragg scattering; compare 
with Ref . 4, pp. 68-69) . 
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The physical justification of the above approximations follows from 
(40) by noticing the role played by the index covariance T in the in- 
tegrand. The weighting introduced by r means that scattering from cen- 
ter y depends upon the neighborhood of y with linear extent L , where 
L is the outer scale. First, the local plane-wave approximation (48) is 
poorest near the origin where the wavefronts are most curved. With a 
criterion of not more than 7r/16 radian departure from plane-wave phase, 
Fig. 2 shows that | y | must be larger than 4Lj;/X. In fact this usual far- 
field condition can be replaced by | y \ > 4L (L /X) i which is less re- 
strictive when L > X. This weaker condition, derived in Appendix B, 
follows from an overbound of the phase error in (48) caused by elimi- 
nating the remainder of (47). Second, the scattering approximation 
(43) is poorest near the observation center ij where the phase (and 
amplitude) of (41) can experience large excursions as r', p range over a 
neighborhood of linear size L . A usual far-field condition is | r\ — y \ > 
4L 2 /X or | ij - y \ > 4(| £ | + L„)7x. Again, when L > X, only a weaker 
condition, 



- y I > 4(U I + l 



o(in^y, 



(53) 



need be met. Condition (53) follows from an overbound of the phase 
error in (43) associated with the remainder in (52), (see Appendix B). 
Strictly speaking, the y-integration in (51) must exclude the near-source 
and near-receiver spheres of radius 4L (L /X) i , and their contributions 
must be evaluated separately. In Section V we give a condition necessary 
for this contribution to be negligible. 




_!_ l q/ 2 

2 



X/32 



yi l /2 

OR 



Fig. 2 — Distance for the plane-wave approximation [in fact, only 4L (L /\) 

required]. 



1/-2 
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Apart from the excluded regions of integration, the validity of ap- 
proximation (51) does not rely upon a "large-scale" condition re- 
quiring the wavelength A to be much smaller than some refractive- 
index scale size. But when such a condition is met, (51) yields both 
a maximum angle of important scattering and a finite volume of 
important scattering. In the approximation (51), the refractive-index 
spectrum (52) is evaluated at the local wave vector, 

k(y) - k.(y). (54) 

Suppose there exists an inner scale l such that for | k \ > 2ir/l the spec- 
trum (52) is negligible. Since the maximum magnitude (54) can attain 
is of order 4ir/\, whereas 2x/\ ^>> 2ir/l , it follows that the integrand of 
(51) is large only for values of y such that 

I KV) ~ UV) I < 2ir/l . (55) 

Under the assumption that | k(y) \ = k = 2ir/\, condition (55) yields 
the maximum angle of important scattering. With \f/(y) the angle 
between k(y) and k,{y), as shown in Fig. 3, we have 

I Ky) ~ k.(y) | 2 = 2kl - 2k\ cos +{y) = 4fc 2 sin 2 ^- (56) 

Then (55) and (56) yield cos f(y) > 1 - X 2 /2l 2 or 21J\ sin $(y)/2 < 1 
or approximately \p(y) < X/l . 

These conditions may be used to find the region of important scat- 
tering. Figure 4 shows cylindrical coordinates with origin at the midpoint 
between transmitter and receiver; there is rotational symmetry around 
the transmitter-receiver axis. With tan \f/ constant, we have 



But 



, , , f - oS tan a + ta n 

tan yp m tan (a + 0) = z r - ^" ( 57 ) 

1 — tan a tan /3 x / 



tan a = -z (58) 

tan 8 = L b ■ (59) 

2~ a 
Algebraic manipulations which include completing a square yield 

a2 + i b + d^i) = idrif- ceo) 
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k-k 




Fig. 3 — Angle of scattering. 

Equation (60) is a circle in the a-b plane passing through the trans- 
mitter and receiver-center locations, Fig. 5. The slope of the curve is 



db 

da 

Thus, near the transmitter, 



b + 



L 



b db 



!+• 



da 



2 tan }f/ 

= tan \f/, 



(-L/2.0) 



and near the receiver 



L 
2~ U 



db 
da 



= tan \{/. 



(L/2.0) 

Also, at the midpoint a = 0, (60) yields 

1 



, = — tan — ■ 
2 Vsin ^ tan fj 2 2 



»-#' x 



Fig. 4 — Coordinates for the region of important scattering. 



(61) 



(62) 



(63) 



(64) 
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L TAN* 



Fig. 5 — Region of important scattering (spherical waves) . 
The condition for important scattering is 

COS ^ > 1 — ^Tj2 — cos ^ . 



(65) 



The volume specified by (65) is enclosed within the surface generated 
by rotating the arc of the circle (60) , with if/ = *, around the trans- 
mitter-receiver axis. This volume lies within the volume common to 
two cones with apexes at transmitter and receiver, each cone with 
half-angle *. 

For the large-scale case, X <<C l ^ L , a condition necessary for (51) 
to be an accurate approximation of the covariance (40) is now apparent. 
The volume of important scattering shown in Fig. 5 must be much larger 
than the volumes in which the integrand of (51) is a poor approximation 
of the integrand of (40). These comprise a near-source cone of axial 
length 4Z^ /2 /X 1 and a near-receiver cone of axial length 4(L -f | £ |) 3/2 /X*. 
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An equivalent condition is seen to be that the transmitter-receiver 
distance must be much larger than the axial lengths of these cones, that 
is to say, 

L= M»4(L,+ U|) 8/ YX*. (66) 

It remains to observe that the covariance expression (40) is itself 
an accurate relation provided that the average field is not severely- 
attenuated by virtue of scattering into the fluctuation field. The at- 
tenuation exhibited in (50) , as evaluated throughout the above region 
of important scattering, must be small; that is to say, 

(Imfc) | 7j |«1, (67) 

where Im k is given by (30). Combining conditions (66) and (67) 
yields an interval for validity of (51). When | £ | = 0, this interval is 

tf'VX* « |u | « (Im ft)" 1 . (68) 

In other words, the transmitter-to-receiver distance must be (i) suf- 
ficiently large so that far-field approximations of the covariance are 
valid, and (w) sufficiently small so that single-scatter perturbation 
approximations are valid. 

V. REDUCTION OF THE INTEGRATION FOR THE COVARIANCE- 
SPHERICAL AND PLANAR CASES 

5.1 Spherical-Wave Case 

The central result of the previous section is the approximation (51) 
of the covariance. The problem remains to evaluate the integral 
specified by this approximation. In this section, we introduce a set 
of coordinates which simplifies the integration, the result being (77). 
Although Section V indicates the extent of the important region of 
integration, the result (77) is equivalent to integration over all space 
rather than over only the important region. Under the large-scale 
approximation, the formula (77) is specialized to (83) and to (85) 
for transverse and longitudinal receiver separations. 

For simplicity, we first observe that (51) can be replaced by an 
expression employing the unperturbed field v rather than the average 
field u c . We need only observe from (7) that 

u e = L~ l j + 0(e 2 ) 

where (Lj) = 0. It immediately follows that (8) can be replaced by 
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where u = Lr r j is the field that would exist in the nonrandom medium 
(c = 0) . In our special case, 

u M _ ex P E fc . 1 r I] 
Uo{r) ~ 4tt | r | ' 

and accordingly 

w .( r ) = -2ek 2 J G(r, r')n(r')u {r') dr' + 0(e 2 ) 
can replace (34). Equivalently, (51) can be approximated by 

uln + ^m*/« _ i 



2/ V 2 

~ AeX I rf?/ (4r)Ty [ r'h-yr ' S(l M ~ fc,(?/) D ' (69) 

where the spectrum S is defined by (52) and where if = y/\y\- That is 
to say, the replacement of u e by u corresponds to the replacement of 
Hv) by k y. 

The volume integration can be carried out with spherical coordi- 
nates which have the receiving center rj as their origin. In such co- 
ordinates, the differential of volume of dfldR R 2 , where R = | tj — y \ 
and d£l is the differential of the solid angle. Since k„(y) is a function 
only of the direction of an element rfn relative to the origin at rj, it 
follows that (69) equals 



4e 2 fc 



f / dQ exp [tfc.-a j dR\\y |" 9 S(\ k D y - k.(y) |) } (70) 



(47T) 4 

where the factor in braces is to be evaluated as a function of R with 
kg fixed. 

The angular integration in (70) will use the coordinates in Fig. 6, 
where 6 = corresponds to the direction of the transmitter. The 
radial integration in (70) will employ the angle \p shown in Fig. 7. The 
argument of the spectrum is the square root of 

I M - k, | 2 = fc 2 (2 - 2 cos 4>) = 4fc 2 sin 2 1- (71) 

The law of sines is 

\y | L R 



sin 6 sin (ir — if) sin (^ — 6) ' 



(72) 
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-*■ | (transverse) 



£ (LONGITUDINAL) 



and thus 



Fig. 6 — Spherical polar coordinates at the receiver. 



, T sin 6 
y \ = L — 



sin yp 

dR = L sin 8 
d\{/ sin 2 \p 



(73) 
(74) 



The radial integral in (70) is thus 

(L sin 9)- 1 f d+ S\2k sin |j , (75) 

and (70) becomes 

j&$® /.' d < f rf * exp <* s /." ^ fi ( 2fc - sin !)■ (76) 

The 6-yp integration is over the triangle {0 ^ ^ *r, ^ ^ ^ ir} - 
{O^^^tt, 0^5^^},so that interchange of the order of integration 
yields 

, 9 wwv * f ^ ^ sin I) f d " f ^ exp ( *'' e ' (7?) 

(Z7r) (47Tly) Jo V «/ •'o Jo 




Fig. 7 — Radial variable i? related to angle if/. 
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Further specialization of (77) is made to the cases in which the 
receiving displacement £ is transverse and is longitudinal, Fig. 6. In 
the transverse case, k s -£ = —k \ $ | sin sin <p, and (77) becomes 

^Bvf f d * *(*• sin !) I! de *•»• i * i sin *> ■ (78) 

In the longitudinal case, k s -£ = k \ £ | cos 6, and (77) becomes 



£^J '*+sfa.Bm*)l*dOe*p{ik. 



£ | cos 0). (79) 



Expressions (77) to (79) correspond to integration over all space, 
rather than over only the region of important scattering. Further ap- 
proximations rely upon the cutoff provided by S(k) for k > 2t/1 , 
where l is the inner scale size and l y> X. For the transverse case, (78) 
becomes 

and in the longitudinal case, (79) yields 

2^2 exp (ik„ U I) / o dx S(x) j o dK exp (-» ^-)- (81) 
The K-integral in (80) can be evaluated in closed form, namely 

f o dK J (k U I) = *J& U I) + 1 *U* I « |)H„(x | g |) 

- \ xJ (x I £ DfftCs | 5 |) (82) 

where H r are the Struve functions. Thus, for the transverse case, (80) 
is 

^g^ [ dx xS(x)[j (x | « |) + | J l( x | ? |)ff.(* | { |) 

- | / (*|£ !)#,(* |*|)J. (83) 
The /c-integral in (81) is related to the Fresnel integrals, namely 

I dK exp \- 1 2kT) = \]T\) w I dt exp { ~ %t ] 



m[*m 



" lS{ ~2k7 



(84) 
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Thus, for the longitudinal case, (81) is 



27r(47rL) 



2 exp 



[*m 



- iS± 



2k t 



(85) 



5.2 Plane-Wave Case 

For the spherical-wave case, the spatial covariance(w,(7j + £/2)«*(i|— 
£/2) ) is given by (77) and its transverse and longitudinal specializations 
(80) and (81) or (83) and (85). By way of contrast, we derive the cor- 
responding expressions for the plane-wave case. 

Approximations (45) and (48) show that the covariance (40) is 
approximately 

(4 + 1)4-|)> 

« 4«V. f dy .fSP'^fl S(\ k(y) - k.(y) |) (86) 
j (47r) i -q — y | 

where k(y) is a constant wavevector, | k(y) \ = k in keeping with 
the interchange of u e and u„, and | u (y) | 2 = 1. Here the integral is 
over the volume of a half-space with k perpendicular to the face. 

With spherical coordinates centered at the receiving center 77, (86) 
becomes 

^h [ dSl exp (ik.'QSQ k-k.\) f dR (87) 

{Zt) J 

where the radial integral has /c,-dependent integration limits correspond- 
ing to the half-space interface. Under the large-scale approximation, 
\/l„ « 1, the radial integral is approximately L. For transverse separa- 
tion, (87) is 

*t r* r 2T ( d\ 

*£ / / d0 <Z? sin 6 exp (-ik \ £ | sin d sin <p)S\2k sin -J (88) 



(2t) 



or 



f dd sin 9J o (k | £ | sin d)s( 2k sin ^ ) ■ 



2. ,.— — ,..— ,-v--,, (89) 

Under the large-scale approximation, with small angles yielding the 
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significant part of (89) the covariance becomes 

~^- rdxxS(x)J {x\S\). (90) 

This expression was obtained by Tatarski [Ref. 4, equation (7.64) ] 
to be equal to twice the correlation function for either the log-ampli- 
tude or the phase fluctuation of the total field. But 

S( x ) = — f dr rY{r) sin xr, (91) 

because (91) is a function of | k j, and 

C* {(r 2 — I t l 2 r I/2 r 2 > I P I 2 

/ dxsin (xr)J B (x | $ |) = l { ' ' U ' ' * ' ' (92) 

J ° lo , r 2 < |f | 2 . 

Substituting (91) and (92) into (90) and changing the variable of 
integration shows that 



(11,(1, + |)u*(r, - I)) « 2«t Jf ^ T(r 2 + I £ 



iM 



(93) 



for transverse separation. This is a central result of much of the 
literature (for example, Ref. 3) ; we have obtained this result in a 
simple and novel way. 

For the case of longitudinal separation, (87) is 

4r itr <>>f / a\ 

*S J I dd dip sin exp (ft. | £ | cos 0)£(2/c o sin - j (94) 



(2*0 



or 



^ | T rf0 2 sin I cos I exp [ft. | £ | (l - 2 sin 2 1) S\2k sin |) • 

(95) 

Under the large-scale approximation, the upper limit of the variable of 
integration k = 2k sin 0/2 can be replaced by infinity. Thus, 



M.l V + 2/ W *V ~ 2 



« ^^ exp (ft, I £ |) J" dx xS(x) exp (-* "4^) ' ( 96 ) 
It does not appear possible to simplify (96) by using (91) together 
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with the sine-transform of the exponential in (96). However, the 
special case 

T(r) = exp (-172© W 

is of interest. Then 

S(x, = (2x) 3/2 / 3 exp(-^72). (98) 

Inserting (98) into (96) and using the variable u = l 2 x 2 /2 for integration 

yields 



(u(v 



+ |)w*(r? - I)) « XLWrf exp (ft. I 6 






(99) 

This expression corresponds to a result of Chernov [Ref. 1, p. 94, equa- 
tion (187)] for longitudinal log-amplitude or phase fluctuations. The 
magnitude of the last factor in (99) is reduced by 5 4 when | £ | = 2k l 2 . 

5.3 Comparison of the Spherical and Planar Cases 

Notice that the relative variance (zero receiver separation) is the 
same for the spherical and planar cases. That is to say, expressions 
(80) and (81) yield the same variance, relative to the spherical-wave 
power (4ttL)- 2 , as do expressions (90) and (96), relative to the unity 
plane- wave power. 

For transverse receiver separation, the planar-case result (90) can 
be compared with the spherical-case result (80) . The weighting of the 
spectral function xS(x) is j (x | £ |) in the planar case; this | £ | -func- 
tion has its first zero at | £ | = 2.4/s with subsequent zeros spaced 
3.1/.T apart. In the spherical case, the weighting is 



- 1 [ oIkJX 

Jn 



> I 



this I £ I -function is a mixture of functions with "periodicity" larger 
than that of J {x \ £ |). Presumably, correlation lengths would usually 
be larger in the spherical case than in the planar case. 

For longitudinal receiver separation, the planar-case result (96) can 
be compared with the spherical-case result (81). The weighting of xS(x) 



is 



cxpl ~* k 
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in the planar case; this | £ |-function has period 4:irk /x 3 . In the spherical 
case, the weighting is 



"/. * iexp (-' 1 2*r') : 



this | £ | -function is a mixture of longer-period functions, and again 
correlation lengths would presumably be larger in the spherical case 
than in the planar case. 

Physical reasoning also suggests that correlation lengths are larger 
in the spherical case than in the planar case. First, compare the regions 
of important scattering. For the spherical case, this region is sketched 
in Fig. 5, where the angle V is given by (65) . For the planar case, this 
region is a cone with half-angle * and axial length L (the transmit- 
ter-receiver separation being replaced by the distance the receiver is 
imbedded into a half-space of random refractive index) , Fig. 8. Com- 
parison of the two regions suggests the the fluctuation field in the 
spherical case is more directive than the fluctuation field in the planar 
case. 

Second, consider the implication of a more directive fluctuation 
field. The directionality function N can be defined by 



Ui 



v + Mv - I 



= J dOexp(tfc.-DJV(fc.). 



A wave in direction k a contributes a correlation exp(ik a -$), and the 
total correlation is a weighted average of such constituents. The form 
(100) is exhibited by (70) in the spherical case and by (87) in the 
planar case. An idealized directionality function would be constant 
with kg -rj above a threshold and would be zero elsewhere. That is to 
say, (100) would be 



2L* 




Fig. 8 — Region of important scattering (plane waves). 



POINT SOURCE PROPAGATION 1153 

~f dtlexpiik.-itj (101) 

where Aft is a small cap on the unit sphere of size 2t(1 — cos 0) ~ 7r0 . 
When £ is transverse, (101) becomes 

[2ir(l — cos 0)] l I I dd d<p sin 8 exp (— ik a | £ | sin sin <p) 

(102) 

(1 - cos©)" 1 f dds'm dJ (k B \ £ | sin d). 

Jo 

Under the small-angle approximation, (102) yields 

2J 1 (k | £ | 0) 



K I £ I 



(103) 



The correlation function (103) is unity at | £ | = 0, is 0.88 at | £ | = 
X/2ir0 ^ 0.16 X/0, and is zero at | £ | ^ 0.61 X/0. 
When £ is longitudinal, (101) becomes 

/* e r 2 * 
[2tt(1 - cos 0)]" 1 / / dd d<p sin exp (ifc„ | £ | cos 0) 

Jo Jo 

= exp (ik. | J |) ' ~ Tufn ' (1 " n C0S 6)] ( 104 > 

K N ' " «fc I £ | (1 — COS 0; 

ri. i , h 1 - exp [-t2fe. 1 £ 1 Aft/4T] 
= GXP <**• I * l } t2*. | £ | Aft/4. (105) 

The correlation exp (ik | £ |) associated with a plane wave is modulated 
by a function having ripple in its numerator with period X(27r/Aft). 

The X-dependence exhibited in (103) and (107) must be tempered by 
the X-dependence of the angle 0. Figures 5 and 8 suggest that would 
be at most \/l (Aft at most ir\ 2 /l 2 ) for equivalence of the idealized and 
true directionality functions. For the transverse case, the null of (103) 
would be at | £ | =0.61 l„ or more; the X-dependence disappears as in 
(80) and (90). For the longitudinal case, the period in (105) would be 
at least 2l 2 J\\ this period is to be compared with the width of the last 
factor in (99) for the plane-wave gaussian-index correlation case. 

Both (103) and (105) exhibit the fact that correlation lengths are 
inversely proportional to Aft, the width of the directionality function. 
But the scattering volumes depicted in Figs. 5 and 8 show that this 
width is smaller in the spherical case than in the planar case. This 
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physical reasoning corroborates the previous interpretation of the 
integration formulas which showed larger correlation lengths for the 
spherical case. 

VI. EXAMPLES OF TRANSVERSE COVARIANCES 

It has been shown that, for transverse receiver separation, the covari- 
ance (u^v + £/2)u*.(ij — f/2)> is given by (80) in the spherical case and 
by (90) in the planar case. These expressions are now evaluated in closed 
form for two illustrative spectra. 

Recall that the spectrum S(k) is related to the refractive index 
covariance r(r) by (52) which becomes (91) for the statistically 
isotropic case. For convenience, (91) is repeated here, together with 
its inverse : 

S(k) = —[ dr rT(r) sin kt (106) 

T(r) = A- f (IkkS(k) sin kt. (107) 

Zir T Jo 

Also, the planar-case covariance (90) is 



^ f dKKS(K)J (K\^\), (108) 

Zir Jo 
and the spherical-case covariance (80) is 

C^ />/.(« !*!)./>«*); (io9) 

the order of integration has been changed. 

The normalization of the spectrum follows from (107) evaluated 
at r = 0, 

1 =±T 2 f"dKK 2 S(K), (110) 

so that c 2 plays the role of the variance of the refractive index. 
Our first example is the case of an exponential spectrum: 

S(k) =7r 2 A 3 exp[-A K ], (111) 

r( f ) = [1 + (r/A)T 2 (112) 

where A ^> X. Then, the planar covariance (108) is 

^-^ [1 + (U I/A) 2 ]" 372 (H3) 
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and the spherical covariance (109) is 



eXLv'A 



[1 + (| | l/A) 2 ]- 
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(114) 



(47rL) 2 27r 

The respective correlation functions in (112) to (114) are plotted in 
Fig. 9. The correlation length for the spherical case is larger than 
the comparable correlation lengths for the planar case and for the 
refractive index. 

Our second example is the case of a constant spectrum: 



SCO = 



r(r) = 



4tt 


t 


^ 


k ^ 2tt/ 


0, 




K > 


2tt/A 


sin 


2irr 
A 


2irr 
A 


2irr 

cos — — 
A 




1 


(2wr\ 


3 




3 


\a) 





(115) 



(116) 







. A = "index correlation length" 
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Fig. 9 — Correlations for the exponential spectrum. 
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Then, the planar covariance (108) becomes 

LMii/aJ (U7) 

The correlation function in brackets agrees with (103) with © = A/A, 
which determines the angular extent of the constant directionality- 
function. The spherical covariance (109) is 

e 2 k 2 L 3ttA \ r fa\Z \ \ _ l \ A 
(4ttL) 3 2 l °\ A / 2x U l/A 

+ .ji^mHi^m - ^(4UK(^, r , (118) 



where H„ are Struve functions. The correlation functions in (116) to 
(118) are plotted in Fig. 10. As before, the correlation length for the 
spherical case is larger than the comparable correlation lengths for 
the planar case and for the refractive index. 

VII. SUMMARY 

In Section II, the perturbation theory of Keller is developed in a 
novel way. 2 This development shows that the nearly uncoupled equa- 
tions (12) and (13) arise naturally from the fundamental pair (3) and 
(4). In Section III, the equation for the average field (12) is solved for 
the case in which the refractive index is statistically homogeneous and 
isotropic. The spherical wave (16) is shown to be a good local approxima- 
tion of the average field; the wavenumber k is a weak function of position 
and satisfies (21). Beyond a few correlation lengths from the source, 
the wavenumber is a constant given approximately by (28) to (30). 

In Section IV, the equation for the fluctuation field (13) is shown to 
imply (40) for the spatial covariance of the field. A useful approxima- 
tion of the covariance is (51) which is then justified on physical grounds. 
For the large-scale case X « l„ , where X is the wavelength and l„ is the 
inner scale for the refractive index, this approximation shows that the 
region of important scattering is given by (60) with (65) and lies within 
the volume common to two cones with apexes at transmitter and re- 
ceiver, each cone with half-angle approximately \/l . The interval for 
validity of (51) is given by (68) which states that the transmitter-to- 
receiver distance is sufficiently large for far-field covariance approxima- 



POINT SOURCE PROPAGATION 



1157 




12 3 4 

NORMALIZED RECEIVER SEPARATION 



!»|l |/ A 



Fig. 10 — Correlations for the constant spectrum. 

tions to be valid but is sufficiently small for single scatter perturbation 
approximations to be valid. In Appendix B, it is shown that a far-field 
condition relative to the covariance is less restrictive than a far-field 
condition relative to the field itself. 

In Section V, the volume integration of (51) for the covariance is 
transformed to the angular integrations exhibited in (77). For the 
large-scale case, the covariance is given by (80) and (81) for trans- 
verse and longitudinal separation of receivers. These expressions for 
our spherical-wave model are contrasted with (90) and (96) for the 
plane-wave model, showing that relative variances are the same but 
that correlation lengths are larger in the spherical-wave case than in 
the plane-wave case. This is to be expected on physical grounds, for 
comparison of the volumes of important scattering for the two cases 
indicates that the fluctuation field is more directive in the spherical 
case. But a more directive field has longer correlation lengths; this is 
illustrated by (103) and (105) for transverse and longitudinal sepa- 
rations under the idealized directionality function in (101). In Sec- 
tion VI, two special cases of refractive-index correlation which cor- 
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respond to an exponential spectrum and a constant spectrum are 
considered. For transverse separation, the covariance functions are 
derived in closed form. Plots of the correlation functions show that 
correlation lengths for the spherical wave case are larger than the 
plane-wave correlation lengths, which are comparable to the correla- 
tion lengths of the refractive index. 

VIII. ACKNOWLEDGMENT 

Discussions with B. H. Bharucha, T. L. Ho, J. N. Lahti, and D. C. 
Stickler are gratefully acknowledged. 

APPENDIX A 

A Mean-Value Theorem 

We show that any solution of 

[V 2 + k 2 ]u(r) = -5(r) (119) 

satisfies the mean-value relation 



-j^2 J dS u(r + p) 



4ttR 

, . sin k R A D - , 

w W^p~ ' <R < \r\, 

lcK (120) 

, N sin kR , sin k(\ r I — R) „ . ■ , 

U{r) -W~ + 4,1, |tB ■ K>|r| ' 

where the integration is over the surface of the sphere {p: \p\ = R}. 
In particular, when u(r) is of the form (16), then (120) becomes (19). 
Introduce a function tj/ (r) that satisfies 

[V s + fc 2 ]iKr) = -8(r- r ). (121) 

Then (119) and (121) imply that 

V-(^Vi* - wW) - u8(r - ?■„) - \f/S(r). (122) 

For the sphere {?*: | r — r \ = R) with the outward unit normal p = 
(r — r e )/\ r — r„ |, the divergence theorem yields 

] dStyfi-Vu - ufi-Vf) 

= >(0, o<a<|r.|, (123) 

V(r tt ) - #0), R > \r„\. 
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We choose i/» to be a linear combination of 

exp (±ik \r - r \) ^^ 

4ltt \r — r \ 

such that $ is zero on the surface of the sphere and satisfies (121). 
This choice is 

m = - l ?y |r ~ r ; l ~i5 ) - < 125 > 

rv ' 4t I r — r | sui kR 

The radial component of the gradient of (125) evaluated on the sur- 
face of the sphere is 

P-W= — A J R rp - (126) 

4xK sm kR 

Then (125) and (126) in conjunction with (123) yield (120). 

APPENDIX B 

An Improved Far-Field Condition 

The kernel (41) used in the integral (40) , yielding the covariance, is 

<*. *m* ,o = ^'p/'ilJ-T'i' ^ (127) 

With the definitions and inverse relations (36) to (39) , 



£ t 

?" = *7 + ^ P = V — o ' 



(128) 



(129) 



r' « y + 2 > p' = U ~ 

the kernel is 



x 



exp [tfc.fl y - y + £(£ - a-)l - U - y - KS - s)l)] n om 
(4tt) 2 |, - j, + «t - ^ I \-n-V-W-x)\ K * } 

The far-field (Fraunhofer) approximation arises from the series 
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expansions (42) 



2 [i? - y 



\v-y-te-x)\ = \ v -y\ - Ir: — K-r-ft - *) + 



(131) 



2 | t, - y 



and the approximation kernel (44) is 



expifcJ-r 2 ^-ft-a) 

L| 17 — y J 



(4x) 2 | t, - y | 2 
A usual condition for the validity of a far-field approximation is 



(132) 



I v - y I » ~ i < 133 > 

where L is an outer scale of the scatering medium. This condition is 
relative to approximation of the field. But relative to approximation 
of the covariance of the field, the condition of validity is 

\v-y\»L [^J. (134) 

In the case L » X, condition (134) is considerably less restrictive than 
condition (133). The reason for this improved state of affairs is that, 
rather than approximating Green's function G, we are approximating 
the kernal GG*. In the computation of the phase of this kernel with 
expansion (131), there is cancellation of terms that ordinarily remain 
when computing the phase of G itself. Overbounding the effect of all 
neglected terms, not just the first one, leads to condition (134). 

Our task is to approximate the phase in (130), namely, the argu- 
ment of the exponential. We put 

Y = v - y, X = S-x, (135) 

so that 

|r-r'| = IF + PM, 

\P- P'\ = \Y -hX\- 
Then, we observe 



(136) 



r±\ I = l'| U±^# + ^y. (137) 
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Next we put 



so that 



Y-X 



a = 



= 



4F 2 



| r - r' | - | Y | (1 + a. + 0)*, 

| p - p' | - | Y | (1 - a + /3)». 



(138) 



(139) 



The next step is to assume | ± a + /3 \ < 1 and expand | r-r* | and 
p—p | with a binomial series. Then, 

(1 ± a + 0)* = 1 + i(±« + P) ~ j£ (±« + i3) 2 

+ ^ (±« + « 3 - i^TT^ (±a + tf + • • • . (140) 



2-4-6 



2-4-6-8 



In the above expression, only terms with differing signs contributed to 
the difference | r—r f \ — | P —p |. Thus, 

(1 + a + /3)» - (1 - a + /3)* 

= a - 2 |^ 2a/3 + 2 |£| « 3 + 2 |^| Zaf + fl, (141) 

where the remainder R has the series expansion 

11-3-5 



R = - 



2-4-6-8 

1-1-3-5-7 



[(a + /3) 4 - (-a + /3) 4 ] 



+ 



[(a + /3) 5 - (-a + P) s ] - 



2-4-6-8-10 

The series is readily "majorized," with the result 

1-1-3-5 



(142) 



R I < 2 



2-4-6-8 



(I « I + | /3 |) 4 



h i'M?'l!'L (|«l + l/3|) 5 + 



2-4-6-8-10 
or, with y = | a I 4- | | and 0< y < 1, 



2 128 



5 4 [\ , 7 
"28 7 L 1 + 1( 



■ 7-9 2 

10 7 i ~ 10-12 7 " i ~ 



(143) 



(144) 
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But the series in brackets is overbounded by 
1 



1 - 7 



= 1 + 7 + 7 + 



(145) 



Thus, y < ^ implies 



^i^i^i + i^ 4 - 



(146) 



The above calculations show that 



r -r' \- \ p - p' 



r 



+ 



I I * I |_ Y a 8 Y 2 

I(¥) 3 +iI¥(fJ +k ]' ci«) 



where 



The conditions 



and 



R l<32 



Y-X 



± 



Y* 

Y-X X 2 
Y 2 + 47 2 



+ 



47V ' 



< 1 



Y-X 



+ 47 2 < 2 



(148) 
(149) 
(150) 



have been imposed. Since condition (150) implies (149) it is clear 
that both are met when 



r 3 



(151) 



The far-field approximation of the kernel employs the leading term 
of (147) ; that is to say, the Fraunhofer phase is 



k -¥—.Y - 2t 3 ~ y .ft - r -) 
k '\Y\' X ~ A |,-y| tt - T) - 



(152) 



The phase error is then 

' X I Y-X X 2 



-k„ 



8 I Y I I X I F 2 



■MnmW-nfi-m*} 



(153) 
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Our task now is to overbound the magnitude of the phase error. 
We impose | X | ^ L„ , an outer scale of the refractive-index correlation 
function I\ This is appropriate for zero receiver separation, | £ | = 0; 
later, L a could be replaced by | £ | + L for nonzero receiver separation. 
Apart from the remainder, (153) is seen to be a function of the cosine 
u-v, with u = X/\ X | and v = Y/\ Y |, and of the ratio w = \ X \/\ Y \. 
We overbound the product in (153) by a product of overbounds in which 
u • v has distinct values, and overbound the second factor by unity (unity 
is greater than both the largest positive value 1 — rsw 2 and the largest 
magnitude, tsw 2 , of negative values, since w < 1). It follows that the 
magnitude of the phase error is overbounded by 

fc„(g^+ | Y\ |fl|). (154) 

We now overbound our previous estimate of \R\. We assume | Y \ 
> L , even though tighter overbounds can be obtained under | Y \ > 
3L as previously assumed. Then, (148) yields 

i«i<lfe + iFy<l(ifn)' (155) 



or 



We use 



^<^(l)W <o - 38 w ] 



*K§(nrrJ- 1)7 



The above phase-error bound (154) is less than 



*-#( 1 + fn^ 



which is less then 



£tf 21 < 5kJA. (m) 

$Y 2 5 < 8 Y 2 (159J 

All of the above calculations required \Y\ > L . With the stronger 
condition | Y | > 4L , the next -to-last calculation becomes 

. Ll ( 16 L„ \ kJZ9 nrm 

k 'Ef\ 1 + TT7\) < sf a 6 (160) 
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This overbound is less than 



4F 2 



(161) 



Suppose we impose the condition that the phase error be less than 
tt/32 and we ask what value of | Y |, the scattering range, is required. 
The last overbound yields 



h - y | = | Y | > 4(^) L e . 



(162) 



This is a less restrictive condition than one specifying the far-field range 
to be much greater than L 2 J\. When L a » X, it is a considerably weaker 
condition. On the other hand, with L > X (but now L ^ X), the con- 
dition (162) still implies the assumption | Y \ < 4L„ under which it was 
obtained. When L < X, the results are valid but vacuous, since the 
pertinent condition is then | Y \ > 4L„ . 

We turn to the approximation of the coherent-field function (46), 



u [y-h 



■•l»-o'- 



exp 



ilk 



y + Z 



- k 



V - o 



(47T) 1 



v + 



y - 



(163) 



in which k is a weak function of position. With the identification y = 
Y, x = X, the previous analysis is applicable. The approximation (48) 
to (50) is 



exp (—2 Im k \ y |) exp 



with a phase error less than 



(4tt) 2 I y 



RekL: 

4 I 2/ I s 



t'(Re k) 



[••] 



(164) 



(165) 



In view of the interchangeability of k and k„ the phase error (165) 
is comparable to 

iW (166) 

Then, a phase error in (164) less than -n-/32 requires that the source- 
to-scattering distance I y \ satisfy 



I y I > 4L (L /X) i . 



(167) 
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Our calculation is similar to one by Lahti and Ishimaru, but the 
calculation (and result) is simpler and the final conditions are less 
restrictive. 10 Simplicity is obtained in part because we use variables 
with symmetric form, (128) and (129), and we overbound a quartic 
remainder rather than modify a cubic remainder. Our quartic-re- 
mainder overbound also yields less restrictive conditions, say | Y \ > 
4L (L /A) i implying a phase error less than tt/32 in comparison with 
| Y | > 7L (L /A)* yielding an error less than tt/10. 
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